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The evolution of the properties of a finite density electronic system as the electron-phonon coupling 
is increased are investigated in the Holstein model using the Dynamical Mean-Field Theory (DMFT). 

We compare the spinless fermion case, in which only isolated polarons can be formed, with 
the spinful model in which the polarons can bind and form bipolarons. In the latter case, the 
bipolaronic binding occurs through a metal-insulator transition. In the adiabatic regime in which 
the phonon energy is small with respect to the electron hopping we compare numerically exact 
DMFT results with an analytical scheme inspired by the Born-Oppenheimer procedure. Within 
the latter approach, a truncation of the phononic Hilbert space leads to a mapping of the original 
model onto an Anderson spin-fermion model. In the anti-adiabatic regime (where the phonon energy 
exceeds the electronic scales) the standard treatment based on Lang-Firsov canonical transformation 
allows to map the original model on to an attractive Hubbard model in the spinful case. The separate 
analysis of the two regimes supports the numerical evidence that polaron formation is not necessarily 
associated to a metal-insulator transition, which is instead due to pairing between the carriers. At 
the polaron crossover the Born-Oppenheimer approximation is shown to break down due to the 
entanglement of the electron-phonon state. 



PACS numbers: 71.38.-k, 71.30.+h, 71.38.Ht, 71.10.Fd 
I. INTRODUCTION 



Electron-phonon (e-ph) interaction play an important 
role in virtually all the materials of present interest. To 
be concrete, e-ph coupling is most likely the driving force 
for superconductivity in magnesium diboridei, in the 
alkali-doped fullerides 2 " and in recently studied interca- 
lated graphite compounds^. Despite the central role of 
electron-electron correlation, also the high-T c cuprates 
are now believed to display remarkable e-ph features- 
witnessed by isotope effects^ as well as spectral and trans- 
port properties 6 . Jahn- Teller e-ph interaction is one of 
the key interactions in the colossal magnetoresistance 
manganites£&, and it may be important in transition- 
metal oxides^. Last, but not least, many different fam- 
ilies of organic materials are characterized by coupling 
with ionic degrees of freedom, from nanotubes to DNA 10 . 
Given the broad variety of different physical properties 
and origin of the coupling, such a wealth of materials cov- 
ers basically all the different regimes characteristic of the 
e-ph interaction. In manganites and cuprates as well as in 
most of the oxides-based materials, the electronic band- 
widths, even if renormalized by strong e-e interactions, 
are larger than phonon frequencies leading to an adia- 
batic character of e-ph interaction, in which the Born- 
Oppenheimer framework is, at least, the reference frame- 
work. This is not the case of highly-polarizable crystal 



lattices of organic semiconductorsii. In this materials 
molecules are bound by very weak Van der Waals forces 
while intramolecular phonons can have quite large oscil- 
lation frequencies leading to an intrinsic anti-adiabatic 
regime of e-ph interaction, in which the adiabatic princi- 
ple breaks down. 

In system with strong e-ph coupling, the carriers lose 
mobility, eventually acquiring polaronic character. A po- 
laron is a state in which the phonon and electron degrees 
of freedom are strongly entangled, and the presence of an 
electron is associated to a finite lattice distortion, which 
in turn binds the electron leading to the so-called self- 
trapping effect. Polarons also tend to create bound pairs, 
called bipolarons. One of the purposes of our studies is 
to clearly distinguish between polaronic and bipolaronic 
features which are often confused in literature. 

The aim of this work is to provide a thorough analysis 
of the Holstein model at half-filing, comparing spinless 
and spinful fcrmions and discussing in detail the role of 
the adiabatic ratio. The backbone of our presentation 
is a numerically exact solution of the Dynamical Mean- 
Field Theory (DMFT) a quantum version of mean-field 
approaches which does not rely on any small parame- 
ter assumption. DMFT has been already applied to the 
study of polaronic systems using analytical result for the 
impurity solution at low densit y 12 ' 13 . At finite density ex- 
act analytical method are not available and we use Exact 
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Diagonalization (ED) as an impurity solver. Such results 
allow us to identify the virtues and defects of various an- 
alytic approximate schemes. In particular we discuss in 
detail a Born-Oppenheimer (BO) scheme which is based 
on the adiabatic limit, and discuss an antiadiabatic ap- 
proach slightly different from the popular Lang-Firsov 
based approaches to polaronic systems^*. In particular 
the BO scheme which will be presented here is suitable 
to be applied to both the weak and the strong coupling 
regime. 

We organize our presentation as follows: After a brief 
introduction of the model and of its treatment within 
DMFT (Sec. [n}, we anticipate the main results from the 
exact DMFT comparing numerical results for different 
regimes (Sec. Illlj) . Then we discuss in detail the Born- 
Oppenheimer approach and compare it with numerical 
results in the adiabatic regime (Sec. HV(1 . and analogously 
we compare the Lang-Firsov approach with results in the 
antiadiabatic regime in Sec. [Vj. Concluding remarks are 
reported in Sec. IVII In appendix are reported the details 
of the calculations used in Sec. IIVI 



where the phonons are defined only on the impurity site 
0, and they interact with the electrons that jump on that 
site. For the Bethe lattice of half-bandwidth D the self- 
consistency enforcing the DMFT solution is given by 15 



II. THE MODEL IN DMFT SCHEME 

The Holstein Hamiltonian is perhaps the simplest lat- 
tice model to describe e-ph interactions. Tight-binding 
electrons are coupled to dispersionless local vibrational 
modes. The Hamiltonian is 

H = -t ( c i* c 3,<? + h - c -) - 9^2(^,* - \){ai + 4)- 



where Cj lCT (cj ff ) and a, (af ) are, respectively, destruction 
(creation) operators for fermions and for local vibrations 
of frequency loq on site i, ni t „ — c\ a Ci tC the electron den- 
sity per spin, t is the hopping amplitude, g is an electron 
phonon coupling. We always fix the chemical potential to 
the particle-hole symmetric value, which fixes the density 
per spin to n = 1/2. In the spinless case there is no sum 
on a. We choose as parameter of the model the electron- 
phonon coupling constant A = 2g 2 /ljqD where D is the 
half-bandwidth of the electrons, and the adiabatic ratio 
7 = Uo/D. 

In DMFT, the lattice model is mapped onto an im- 
purity problem subject to a self-consistency condition, 
which contains all the information about the original lat- 
tice. Our model Q becomes a Holstein impurity model 
(HIM), 
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where G(iu> n ) is the local Green's function of the system. 
One could eventually get the electron self-energy trough 
the following relation which holds in the Bethe lattice 
case 



G(w) 



££G(u) - E(w) 



(4) 



We solve the HIM by means of exact diagonalization, 
i.e., by truncating the sums in the first two terms of Eq. 
(J2J) to a small number of terms Nj,, so that the Hilbert 
space is small enough to use, e.g., the Lanczos algorithm 
to compute the T = Green's function. For the case 
of phonon degrees of freedom we consider here, also the 
infinite phonon space has to be truncated allowing for 
a maximum number of excited phonons N p h- In all the 
calculations presented here the convergence of both trun- 
cations have been checked. The value of N p h has to be 
chosen with special care in the adiabatic regime and in 
strong coupling, where phonon excitations are energeti- 
cally convenient. As far as the discretization of the bath 
is concerned, the convergence of thermodynamic averages 
and Matsubara frequency properties is exponentially fast 
and Nj, ~ 8—9 is enough to obtain converged results. The 
method also offers the advantage of a direct evaluation 
(1) of real- frequency spectral properties such as the electron 
and phonon spectral functions. The main limitation is 
that these quantities reflect the discrete nature of our 
system much more than their imaginary-frequency coun- 
terparts. In practice, the spectra are formed by collec- 
tions of 5-functions. Of course this limits our frequency 
resolution, and suggests that the method is better suited 
to gain knowledge on the main features of the spectra, 
rather than on the fine details. 

In the following we define some important quantities 
which we use to discuss and characterize the physics of 
our model. Our focus is mainly on the polaron crossover 
and the metal-insulator transition. The quasiparticle 
residue Z = (1 — d'E(u))/du))\~_ , E(u;) being the electron 
self-energy, proved extremely useful in marking metal- 
insulator transitions, both for repulsive and attractive 
models, and it has already been employed for the char- 
acterization of the bipolaron transition^. More detailed 
informations on the passage from a metal to an insulator 
can be gained through the study of the electron spectral 
function (or density of states) 



p(lu) = ImG{ui) 
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As far as the polaron crossover is concerned, a key 
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quantity is the lattice displacements probability distri- 
bution function (PDF)AL±£ 

P(X) = (V> |X)TOo>, ( 6 ) 

where \ipo) is the groundstate vector of the impurity 
model and X is the phonon displacement operator on 
the impurity site. The appearance of a bimodal PDF 
signals indeed at finite densities the formation of a pola- 
ronic state, i.e., a state in which the presence of the elec- 
tron is associated to a definite polarization of the lattice. 
At finite densities and for strong coupling the electrons 
are able to drastically change the phonon properties, as 
opposed to the case of a single polaron 17 . In the spinful 
case such a state may be bipolaronic, i.e., the polariza- 
tion of the lattice may be associated to a local pairing of 
a pair of electrons with opposite spin. 

Further information on the polaronic properties and 
on the mutual interaction between electrons and phonons 
can be extracted from the phonon propagator 

D(t) = -i (T{a\t) + o(t))(ot(0) + o(0))) , (7) 

from which the phonon spectral function is readily de- 
fined as 

p P h(u) = ImD(w), (8) 



III. RESULTS 

In this section we briefly anticipate the main results 
obtained through ED solution of the full DMFT equa- 
tions. We compare spinless and spinful fermions orga- 
nizing the results according to the degree of adiabaticity. 
The results are summarized in Figs. ^and|2 respectively 
referring to the adiabatic and anti-adiabatic regime. In 
each figure the left column are dedicated to the spinless 
fermion case, and the right one to the spinful system. 

As we discuss in Sec. IIVI in strong coupling, we can in- 
troduce relations between the parameters A and 7 which 
allow for a direct comparison. For the adiabatic param- 
eter we have J sp inful = 2 'spinless- We choose therefore 
7 = 0.1 (7 = 0.2) as representative of the spinless (spin- 
ful) adiabatic regime and 7 = 2.0 (7 = 4.0) for the an- 
tiadiabatic regime. 

Let us now qualitatively describe the behavior of the 
model in these two regimes by comparing the phonon 
PDF (First row in Figs. □ 0, phonon DOS (Central 
Row in Figs. ^ 0) and electronic DOS (Bottom row in 
Figs. Let us remark again that the discretization of 

the electronic Hilbert space of the bath allows to extract 
only gross features from the spectral properties of both 
electron and phonon. 

We first discuss phonon properties through an analysis 
of the evolution of the phonon PDF and DOS. The PDF 
has a qualitatively similar behavior in adiabatic and an- 
tiadiabatic cases for spinless and spinful cases: increasing 



e-ph coupling the PDF becomes bimodal signaling a po- 
laronic or bipolaronic phase 16 ' 17 . However the value of 
e-ph coupling for which the bimodal behavior establishes 
is much higher in the antiadiabatic case. 

The phonon DOS shows instead a rather strong qual- 
itative dependence on the adiabatic parameter. Ap- 
proaching the polaronic region in the adiabatic regime 
we observe a softening of the phonon frequency. In the 
anti-adiabatic case, instead, we have a transfer of spec- 
tral weight from high to low frequency while the position 
of the resonances are roughly unchanged. In this regard 
a qualitative difference between the spinless and the spin- 
ful case appears at low energies: in the former case the 
low energy peak disappears at strong coupling. 

The most important qualitative difference between the 
spinless and the spinful case is the behavior of the elec- 
tronic DOS, and it is more evident in the anti-adiabatic 
case. In the spinful case a sharp Metal-Insulator Transi- 
tion (MIT) transition, signaled by the opening of a gap 
in the electron DOS, is observed for a critical value of 
the coupling which is not dramatically dependent on the 
adiabatic parameter. Conversely in the spinless case no 
sharp MIT takes place, even well inside the polaronic 
region. 

All the above observations can be summarized in the 
phase diagrams of Fig. [21 The points in which the 
phonon PDF becomes bimodal are used to draw a line 
which represents an estimate of the polaron crossover re- 
gion. This line is apparently strongly 7 dependent in 
both the spinless and spinful case. Beyond this line a 
polaronic (bipolaronic) regime is attained in the spin- 
less (spinful) case. In the spinful case, we can also 
draw a MIT line, which separates a normal phase from 
a paired insulating phased 6 ". Notice that in the antia- 
diabatic regime, where the Holstein model is approach- 
ing a purely electronic attractive Hubbard model, a pair 
has not necessarily a well definite associated polariza- 
tion therefore at very large value of 7 we can have pairs 
without bipolaronic behavior even in the Holstein model. 
In the following sections we will analyze separately the 
adiabatic and the non adiabatic regimes of the system, 
discussing the numerical results in comparison with dif- 
ferent analytical approaches suitable for the two opposite 
regimes. 



IV. ADIABATIC REGIME 
A. Adiabatic limit and BO approximation 

We briefly describe, as a starting point of the Born- 
Oppenheimer (BO) procedure, the adiabatic limit in 
which 7^0 keeping A fixed. This limit has been thor- 
oughly studied in Reftf, and here we introduce a slightly 
different, yet equivalent, formulation. 

In the adiabatic limit the phonon displacement be- 
comes a classical variable, therefore the HIM |J2J can be 
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FIG. 1: (color online) DMFT data in the adiabatic regime 7 = 0.1 spinless (panels on the left) and 7 = 0.2 spinful (panels 
on the right). In each panel the various Curves refer to different value of A spanning from 0.1 to 1.8 in the spinless case and 
from 0.05 to 1.1 in the spinful case and are shifted according A value. The first line show the phonon PDF, the central line the 
phonon DOS and bottom line the electronic DOS. 

written as oscillator characteristic length. 

. This Hamiltonian represents an impurity electron 

— g'n Q X H — kX 2 . (9) which undergoes multiple scattering, as depicted in Fig. 

^ 0]a), and it can jump on the conduction band via the 

where g' = g£/\/2 and £ — y^h/MuJo is the harmonic first term of @. The Green's function is immediately 




FIG. 2: (color online) DMFT data in the antiadiabatic regime 7 = 2.0 spinless (panels on the left) and 7 = 4.0 spinful (panels 
on the right). In each panel the various curves refer to different value of A spanning from 0.4 to 6.5 in the spinless case and 
from 0.2 to 3.0 in the spinful case and are shifted according A value. The first line show the phonon PDF, the central line the 
phonon DOS and bottom line the electronic DOS. 



written as weights 19 20 . These values can be obtained by minimiz- 

ing the total ground state energy which we call adiabatic 



G(o,) = ^>, (10) 
1 Go ( UJ ) + 9 x i 



where I labels all possible value of X giving the same 
ground state energy and wi are the corresponding 
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FIG. 3: (color online) Phase diagrams of the spinless (left) 
and spinful (right) Holstein model at half filling. Solid lines: 
numerical results from DMFT, dashed lines: approximations. 
Left panel: the bold line is the polaron crossover from bi- 
modality of P(X) and the dotted line is the anti-adiabatic 
estimate a 2 > 1 for the polaron crossover. 
Right panehbold curve is the bipolaronic MIT from vanishing 
of Z, thin solid line the polaron crossover, bold dotted line is 
the anti-adiabatic prediction for bipolaronic MIT (Eq. 1371 1. 
light dotted line is the anti-adiabatic estimate a 2 > 1/4 for 
the polaron crossover. 



potential' 



V(X) = l -kX 2 - 9 -\X\ 



iuj n + g'X 



(11) 



where s is spin degeneracy. 

In formula Eq. (jl 1|> we have found useful to separate 
the contribution in absence of hybridization ( first line of 
Eq. dJ) from a remainder (last line). The latter term 
can be obtained easily trough the linked cluster theorem 
as depicted graphically in fig^]b) 



FIG. 4: a) The equation for Go (thin line) and G (bold line). 
Dashed line is the single site impurity propagator (1/w) bold 
arrow is the hybridization constant 14 double line the bath 
propagator and x-type insertion the scattering with static dis- 
placements field, b) The diagrams expansion of the adiabatic 
potential. 

The equation for the extrema of the adiabatic potential 







g'X„ 



(12) 



The self-consistency condition © together with Eqs. 
(|10|l and i|12|) . completely solves the problem. Notice 
from Eq. H12|) the correspondence between the spinless 
and and the spinful case upon rescaling A = g l2 /kt to 
A/2 in the latter case. 

In the Bethe lattice case at half filling it can be shown 
that for A < A c , where A c = 37r/(8s), Eq. i(T2|l has only 
one solution while for A > lambda c , three solutions exist: 
two stable (X[ ^ and wi = 1/2) and one unstable 
(X = 0) 17 . Eq. ((TU|> can be recast as an average over the 
phonon PDF which is a single ^-function for A < A c and 
splits in two symmetric ^-functions for A > A c . 

A MIT occurs at a larger coupling Xmit = 1.328/s as 
a vanishing of the DOS at the Fermi levelit 

As a final remark we mention the similarity be- 
tween the adiabatic limit of the Holstein model and the 
DMFT solution of the Falicov-Kimball model 22 as al- 
ready pointed out in Ref. 20 . This similarity is evident 
once we consider the Coherent Potential Approximation 
(CPA) form of the Green's function (|10fl and compare it 
with, e.g., that given in Refw 

The BO procedure goes on by quantizing the adiabatic 
potential after adding the phonon kinetic energy contri- 
bution. Introducing the scaled variable u — g' Xj D\fs\ 
the BO phononic hamiltonian reads 



Hbo 



7 d 



sV{u). 



(13) 



and V(u) is given by Eq. (|11|) in the variable u. Notice 
that the spinful BO hamiltonian maps onto twice the 
spinless one upon rescaling A/2 — > A and 27 — > 7. 

Whether the gradient term represents the most rele- 
vant contribution from quantum fluctuation of phonons 
can be questionable. A more general non-local contri- 
bution to the adiabatic potential arises in the effective 
phonon action as it is discussed in Refi^4. As it is dis- 
cussed there this non local part may affects the determi- 
nation of the phonon properties. However we decide to 
pursuit the way of simplicity and discuss the BO approx- 
imation in view of comparison with ED data. 

While phonon properties are immediately obtained at 
this stage from the solution of the one-dimensional an- 
harmonic system of Hamiltonian Eq. i|13|) . electronic 
properties must account non-trivially for the tunneling 
of phonon coordinates. 

The simplest way to describe electrons coupled to a 
tunneling system is to map it onto a two level system. 
In our model this can be accomplished by changing the 
basis (operators a) from that of the harmonic oscillator to 
the that defined by the solution of i|13|) . Then projecting 
out all the states but the first two (|+), |— )) we get the 
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following two state projected model (TSPM): 

(7 ' k,CT 

+ £^(/K*+4, ff /*)> (14) 

where a x and a z are two Pauli matrices in the space 
spanned by |+) and |— ) and the quantities e and A are 
given by 

e = ff £(+|o + o t |-> (15) 
A = ^«+| a t a | + (_)_| t a |_» (16) 

The latter quantity A as a clear meaning as a tunneling 
frequency between the two phononic states. A similar 
model was introduced in RefiSi to study the strong cou- 
pling limit of the Holstein model. Here we remark that 
we use this model as a tool to gain physical insight in the 
analysis of numerically exact output of DMFT. 

A sample plot of the parameters of this model as ob- 
tained from the adiabatic limit of the Hamiltonian Q is 
reported in fig. |3J At weak coupling we obtain e ~ g~ 
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FIG. 5: (color online) Parameters of the TSPM in the spinless 
case for 7 = 0.1. The spinful case is simply obtained by taking 
^spinfui = ^spinless / '2. The adiabatic polaronic transition A c 
is marked by a solid arrow while the adiabatic MIT is marked 
by a dashed arrow. 

and A ~ At strong coupling instead e scales as A/4 
while A vanishes exponentially, even if it never becomes 
strictly zero for any finite value of e-ph coupling in both 
spinless and spinful case. 

There are two limits in which the TSPM reproduces 
exactly the original Holstein model within DMFT: the 
weak coupling and the adiabatic limit. In the former 
case the projection of the phonon space has no relevance 
therefore the TSPM reproduces the perturbation expan- 
sion developed (in the limit of infinite bandwidth) in the 
classical Ref. 123. The adiabatic limit is instead recovered 



as A — > 0. No phonon tunneling occurs and the model 
can be solved exactly by CPA recovering the solution of 
Ref. E3 

1 | 1 
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FIG. 6: (color online) The spectral function of the spinless 
TSPM in the zero hybridization case (bold lines) and in the 
CPA approximation (thin line). The upper panel refers to a 
typical weak-coupling situation (A = 0.1,7 = 0.1) while the 
lower panel presents a strong-coupling case (A = 1.6,7 — 0.1). 
In the inset of the lower panel -/mS is shown. 
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FIG. 7: (color online) The spectral function of the spinful 
TSPM in the zero hybridization case (bold line) and in the 
CPA approximation (thin line). The upper panel shows a 
weak-coupling result (A = 0.2,7 = 0.1), the lower panel a 
strong coupling one (A = 1.9,7 = 0.1). In the inset of the 
lower panel — JmS is shown. 

To study the effect of phonon tunneling events in 
strong coupling it is useful to consider the TSPM in ab- 
sence of hybridization (V4 = 0). Taking into account the 
self-consistency condition ©, this is the atomic (D = 0) 
case for the lattice Hamiltonian 0). Let us first consider 
the simplest spinless case where the eigenstates are given 
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by (see also Appendix for details) 



\v+) = |0) ® 



1 



A 



l«a>-/ + |0)« 



y/2l{l - at) \' 
1 



la 

e + la 
A 



(17) 



where a = ±1 labels the phonon "spin" and |0) is 
the Fock vacuum for the impurity. The two-component 
spinors live in the space defined by |+) and |— ). The 
eigenvalues are in this case doubly degenerate: 



= a0l 



(18) 



where (3, a = ±1. and I 2 



A 2 . 



The eigenstates \v~) and \v + ) describe respectively the 
presence and the absence of an electron on the impurity, 
and they are associated to different phononic states. The 
phononic tunneling A allows for the existence of quan- 
tum defects on the impurity i.e. states in which the im- 
purity is occupied (empty) and the deformation is the 
one adiabatically associated to the absence (presence) of 
the electron, e.g., \vf) and \vZ) ■ The adiabatic case is 
recovered once A = so that only |tTt) and \vT) survive. 

The atomic (Vk = 0) Green's function can be easily 
found to be 



G a {0) 



2 h 



i 



P + 2lc 



A 2 1 

1?0 



(19) 



G a {u>) has a pole at u — induced by phonon tunneling, 
whose weight in fact vanishes as A — > 0, accompanied by 
two resonances at ±21, as depicted in Fig. The zero 
energy peak is due to transitions in which both charge 
and phonon " spin" change while the side peaks take into 
account charge transfer in a frozen phonon " spin" . In this 
sense the side peaks are adiabatic features which survive 
when phonon tunneling A — > 0. 

In the spinful case denoting with \v) the singly occu- 
pied states and with \u) the doubly occupied or empty 
states we obtain (see Appcndi>{XJ) 



|«+) = |0) ® 



IO = / T + /f |o) 



A 



y/2l{l - at) V e ~ la 

1 fe + la 

^21(1 + ae) V A 



(20) 



\4) = fp\Q)® 



V2K- 1 



K) = #|o>® 



V2 V 



(21) 



The eigenvalues are 



E^(v)=a(3A 



(22) 



each level is doubly degenerate and (3, a = ±1. 

Notice that the empty or doubly occupied impurity 
states \u) are associated to a finite deformation, while 
the singly occupied states are associated to a tunneling 
phonon state (eigenstates of a x ). 

In this case a single particle excitation has always non 
zero energy as it can be seen by the Green's function 



G o (0) 



1 



1 



2^2/ 

a=± 



I- A 



l + A 



, + a(l + A) + a(l-A) 



depicted in Fig. [7] The most striking with respect to 
the spinless case (|19fl is the absence of the zero frequency 
pole. In the spinful case the tunneling of the phonon A 
is always associated to a finite energy transition, and it 
only splits the finite frequency poles associated to the 
transition from singly to the empty or doubly occupied 
ones. 

To analytically span from the strong (Vk — > 0) to the 
weak (g — ► 0) coupling regimes of the equivalent impurity 
Hamiltonian it is useful to devise a Born-Oppenheimer 
Coherent Potential Approximation (BOCPA) scheme. 
Starting from the Green's function for Vk — (|19|l for 
the spinless case and <|23[1 for the spinful one, we notice 
that in both cases G a {0) can be written ad sum of two 
contribution G a {0) = (l/2)(G a ,+ (0) + G a ,-(0)) where 
(±) label a phonon state as in (|17fl and (|20p. Then we 
write the propagator in presence of hybridization as 



(24) 



where £i& r = J2k ^fc /( w — Ek) and takes into account 
the hopping of the electron from the impurity to the 
conduction states (see Fig. 0J). BOCPA assumes that 
the tunneling states of the phonon in the adiabatic po- 
tential remain unaltered during an hybridization event. 
The results depicted in figs. and [7| As in standard 
CPA at each local level is associated a band but in con- 
trast with the CPA for the Hubbard model, our BOCPA 
gives a Fermi liquid solution in the spinless case for every 
value of the coupling. The low-energy band arising from 
the zero energy pole in the zero hybridization limit is in- 
deed coherent. This can be easily realized by analysis of 
the self-energy. When u> — > the propagator the spin- 
less propagator defined by l|19f) and l|24|l is dominated 
by the zero-energy pole of G a (jT§|l and consequently the 
self-energy obtained through Q is purely real. 

Conversely in the spinful case a MIT transition similar 
to that of the Hubbard model in the CPA approximation 
is observed at a critical value of In both spinless 

and spinful case the weak coupling limit of BOCPA (thin 
lines of upper panels of Figs. |H|and[7I) bears strong re- 
semblance with the classical result of Ref. Finally 
we emphasize that we recover the adiabatic solution of 
Ref. [It] as 7 — > for finite A as A — > 0. In this case 
the BOCPA is exact and gives the Green's function of 
Eqs. (|10fl . However the BOCPA procedure is certainly 
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affected by serious problems approaching the MIT in the 
spinful case. A more careful treatment of the low-energy 
part of the Green's function has been performed in this 
case in Refi^i, where it has been observed a MIT sce- 
nario similar to the half- filled the Hubbard model, i.e., a 
quasi-particle peak that shrinks to zero width approach- 
ing a critical value of A, as it will be discussed at the 
end of the next section. In the spinless case instead a 
resonance is present at zero energy in the even within a 
CPA approach. It is not associated to a Kondo effect but 
rather to phonon tunneling which drives charge fluctua- 
tions. On the other hand a Kondo like behavior can be 
ascribed to the bipolaron or pair formation (spinful case) 
and should be treated either numerically as in our case 
or with a more appropriate theory as in Ref. 24. 

Now let us briefly discuss which kind of processes are 
left out from a BOCPA theory. For the sake of clarity 
here we only present the main results of a perturbative 
treatment up to second order in the hybridization T4, 
and we refer to Appendix 1X1 for more details. 

) = ) + ) > > ) 

(1) 

+ --)--- > ► -->-- 
(2) 

FIG. 8: The second order perturbation expansion in Vk for the 
impurity propagator. Symbols are defined in the Appendix lAl 

The electron Green's function can be written at second 
order in hybridization Vk as G(uj) = Gi(u) + G2(u) where 
G\ is a term which is resummed in BOCPA scheme and 
G2 is an extra term coming from processes in which the 
phonon "spin" changes during the hybridization process. 
The two terms are represented diagrammatically in Fig. 
GO G 2 reads: 

G 2 (w) = £ffaMX£HffaM (25) 

a 

where 

E» = £e(-aE fc ) (26) 

^ uj - Ek + 2la 

k 

and the step function Q(—Ek) represents the Fermi func- 
tion at zero temperature. In Eq. I|25f) 9a (<*>) is in the 
spinless case 

M = f(— l -^r--) (27) 

r \oj + 2la cu ) 
while in the spinful case is given by 

= 1 ( I 1 I a \ - T 7i at ) • (28) 
I \uj + a(l + A) u + a(l — A) / 

As discussed in^ the g propagators do not carry charge, 
as they are correlation functions of the phonon " spin" , 



and we will refer to them as /— spinons. By inspecting 
the structure of the Dyson's equation we see that 
the insulating character is preserved after the inclusion 
of G2 terms in the spinful case. In the spinless case a 
modification of the zero energy pole arises from Eq. I|25() 
but it is still present also after inclusion of the /— spinon 
terms. We notice however that approaching zero en- 
ergy the /—spinon propagator becomes a independent 
in the spinless case. Assuming particle-hole symmetry 
^ Q E^,(o; = 0) = which leads to a vanishing contri- 
bution of the /—spinon diagrams at zero energy leading 
to the validity of the Luttinger's theorem even within 
BOCPA for the spinless case. 



B. Discussion of the results in the adiabatic regime 




-3-2-1 1 2 3-3-2-1 1 2 3-3-2-1 1 2 3-3-2-1 1 2 3 
X X X X 



FIG. 9: (color online) The phonon PDF in the spinless (upper 
panel) and spinful (lower panel) cases. Solid curves are DMFT 
results and dashed lines the BO approximation. The values 
of A are in the spinless (spinful case) from left to right: A = 
0.5(0.25),A = 1.5(0.75),A = 1.7(0.85), A = 2.2(1.1) 

Fig. ^presents a comparison between the BO approxi- 
mation and the DMFT results for the phonon PDF. The 
anharmonicity due to e-ph interaction increases as the 
coupling increases leading first to a non-gaussian and fi- 
nally to a bimodal PDF. This behavior signals the ap- 
pearance of static distortions, even if we are neglecting 
any ordering between them. We can estimate a "cen- 
tral line" of the region in which the polaron crossover 
takes place according to the values of A and 7 at which 
phonon PDF becomes bimodal. From Fig. [5] is evident 
that BO approximation works well in both the metal- 
lic and the polaronic regimes. The reason for the ac- 
curacy of the BO procedure in the polaronic regime is 
that, contrary to its usual implementation in the weak 
e-ph coupling 2 ^, here we take into account the anahar- 
monicity through Eq. in a non perturbative way. 

However, BO does not accurately reproduce the phonon 
PDF around the polaron crossover. In this case electron 
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and phonon states are strongly entangled, and cannot be 
approximated properly by a disentangled BO state. By 
a comparison of the spinless and spinful cases in fig. El 
we see that the occurrence of the MIT does not influence 
much the differences between full DMFT and BO, which 
are in both cases relevant near the polaron crossover. 
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FIG. 10: (color online) Phonon DOS in the spinless (upper 
panel) and spinful (lower panel) cases. Solid curves are DMFT 
results, bold impulses the BO approximation. The values of 
A are the same as in Fig. 

This is not the case of the phonon DOS shown in Fig. 
ITTI In this case a qualitative difference between the spin- 
less and the spinful case appears at low energy as already 
noticed discussing Figs. Q] an dEl hi the spinful case and 
above the MIT (rightmost picture in the lower row of 
Fig llOfl . the low-energy peak of the phonon DOS com- 
pletely disappears. This behavior cannot be predicted by 
BO approximation which gives undamped phonon peaks 
approaching zero with increasing spectral weight. An 
explanation of this behavior needs an anti-adiabatic ap- 
proach which will be described in the next section. 

The data for the phonon DOS appearing in Figs. 2] 
and 1101 can be qualitatively compared to those obtained 
by Numerical Renormalization Group (NRG) solution of 
the impurity model (See, Fig. 5 of Ref. l45|) . The main 
differences between our ED results and NRG seems to 
be the pole at ujo which is present in NRG for all the 
coupling strengths shown, while in our analysis this pole 
softens and finally disappear at the MIT. We ascribe this 
different behavior to the effect of the high energy part of 
electronic spectrum, which is treated approximately by 
NRG but it may strongly influence the phonon proper- 
ties. This is easily realized by analyzing the simple bub- 
ble diagram for the phonon sclf-cncrgy, in which both 
low- and high-energy electronic scales contribute in the 
same way the low-energy part of the phonon self-energy. 

In Fig. ^2 we compare the electronic DOS from full 
DMFT with BOCPA. The different behavior of the spin- 
less and spinful case is not so evident in this adiabatic 
regime. However, comparing the spinless spectrum for 
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FIG. 11: (color online) DOS in the spinless (upper panel) and 
spinful (lower panel) cases. Solid curves are DMFT results, 
dashed lines the BOCPA approximation. The values of A are 
the same as in Fig. |5] 



A = 1.5 with the corresponding spinful for A = 0.75, a 
quasiparticle peak is present in the former case, while a 
depletion of low energy part is much more evident in the 
latter. At strong coupling the discretization of the bath 
degrees of freedom inherent to the ED solution of DMFT 
does not allow us to identify a well defined quasiparticle 
peak in the spinless case. However a more careful analy- 
sis of the quasi particle spectral weight^ shows that no 
MIT occurs in the spinless case. Notice that the BOCPA 
seems to be much closer to DMFT-ED in the spinless 
than in the spinful case where the CPA approximation 
for the electronic degree of freedom is apparently much 
less adequate. 



V. ANTI-ADIABATIC REGIME 
A. Canonical transformation analysis 

While in adiabatic limit the phonon displacement be- 
comes a classical variable, and we are left with an elec- 
tronic model which depends parametrically on it, in the 
opposite limit (7 >> 1) the roles are exchanged, and 
we have a parametrically fixed electronic charge on a 
given site. In this regime the most reasonable starting 
point is the Lang-Firsov (LF) canonical transformation 14 
U — exp(S), which here we only apply to the equivalent 
impurity Hamiltonian (0). The generator of the trans- 
formation reads 

S=-aY,(tiU-\)(a!-a), (29) 

a 

introducing the parameter a = g/ujQ which is the 
relevant e-ph coupling parameter in the anti-adiabatic 
regime. 12,28 Phonon and impurity electron operators are 
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transformed in the following way: 

5-exp(S)aexp(-S) = a + ]T(/t /(J _ I) (30) 

f a = exp(S)f a exp(-S) = f a exp(a(a^ - a)) (31) 

This canonical transformation diagonalizes the impurity 
Hamiltonian in the absence of hybridization by elmini- 
nating the e-ph interaction part. In the spinful case the 
phonon energy term of @ gives rise to the well known 
bipolaronic instantaneous attraction. The hybridization 
term of © is modified according to lj3*Tl) by acquiring an 
exponential term in the phonon coordinates leading to 



e b He 



E 



a(a T -aW (J 



V k (clJ a + h.c.) + 

V E k c\ c k>a - 2— (s - l)n T n[ - 
, — 

k.a 

— V](- - n CT ) + woa f a, (32) 

U!q — ' 



where n a = flf„. 

Notice that in the anti-adiabatic limit 7 — > 00, if A is 
kept constant a vanishes. In this case spinless electrons 
results unrenormalized while spinful electrons are de- 
scribed by an attractive Hubbard model with \U\/D = A. 

If we want to proceed with analytical methods, the 
hybridization term must be treated in an approximate 
way. Assuming that in the anti-adiabatic limit the im- 
purity density is constant during the fast motion of the 
phonon we average out the phonon term on the displaced 
phonon ground state. This is the so-called Holstein Lang- 
Firsov Approximation (HLFA), which has not to be con- 
fused with the exact canonical transformation treatment 
(|29|l . HLFA gives rise to the exponential renormaliza- 
tion of the hybridization constants where each V k of Eq. 

is replaced by V k exp(-a 2 /2). Through the self- 
consistency condition J3J), such a replacement implies the 
well known exponential renormalization of the bandwidth 
D exp(— a 2 ). 

The main results of HLFA can be summarized as fol- 
lows: 

i) In the spinless case the HLFA on the Holstein im- 
purity Hamiltonian generates a non interacting impurity 
which is connected with the bath through an exponen- 
tially reduced hybridization; 

ii) In the spinful case the impurity site presents an 
attractive instantaneous interaction term of the Hubbard 
type U = -2g 2 /aj . 

To get the electron Green's function G(w), the explicit 
action of LF transformation (|30|l and l|31l) have to be 
taken into account. Following Refs. l29land l3Cl we obtain 
in both spinless and spinful cases 



,2|< 



-G p (u> — nu>o). (33) 



where G p (uj) is the Green's function of an impurity (with 
a negative U interaction in the spinful case) hybridized 
with an exponentially reduced hybridization to a bath of 
conduction electrons. The self-consistency condition can 
be written explicitly in the spinless case due to the lack 
of interaction terms on the impurity: 



G„(w) 



1 



(34) 



where G{0) is the local Green's function of the lattice. 
In the spinful case a Lang-Firsov Coherent Potential Ap- 
proximation (LFCPA) can be devised for resulting HLFA 
negative-U Hubbard model giving 



1 



1 



^-e- a2t iG(uj)-U/2 
1 \ 



4 



G(tu) + U/2 



(35) 



A comparison between the zero hybridization (atomic 
D = 0) DOS and results from HLFA and LFCPA results 
is depicted in figs. El El 
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FIG. 12: (color online) The spectral function in the spinless 
case in the zero hybridization case (bold line) and in the HLFA 
approximation (thin line) . The upper panel refers to a typical 
weak-coupling situation (A = 2.7,7 — 1-2) while the lower 
panel shows results for intermediate coupling (A = 5.4, 7 = 
1.2). 

Notice that the theory developed here for the Holstein 
impurity model differs significantly from that developed 
directly in the lattice models. In that case an equation 
identical to (|33ll is recovered for a band of free electrons 
therefore giving a low-energy coherent polaronic band in 
the spinless case. It is however easy to show that this 
form of the spectral function is not compatible with a 
/c-independent self-energy. Self consistency condition eq. 
(|34|l which enforces a local self-energy give rise to a non 
zero damping at the Fermi level even in the spinless case. 
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FIG. 13: (color online) The spectral function of the spin- 
ful case in the zero hybridization case (bold line) and in the 
LFCPA approximation (thin line). The upper panel is dedi- 
cated to weak coupling (A = 2.7, j — 1.2) and the lower panel 
to intermediate coupling (A = 5.4, 7 = 1.2). 



However as far as 7 becomes larger a gets smaller repro- 
ducing the anti-adiabatic coherent behavior at low energy 
in the spinless and the negative-[/ behavior in the spinful 
cases. 

In the anti-adiabatic regime the HLFA approach gives 
an estimate of the MIT by the simple scaling 



Am/t = \U/D\ M iTexp(-a 2 ) 



(36) 



where \U/D\mit = 2.94 is the MIT value of the negative- 
Li Hubbard model 31 which is equal to that of the repulsive 
mode! 32 . As a function of the adiabatic ratio we get: 



1MIT = ~\U/D\ 



A 



MIT 



2\og{X/\U/D\ M iT) 



(37) 



which indicates that the MIT takes place at lower values 
of A as the adiabatic regime is approached (see Fig. 
Of course the previous formula has not to be applied in 
the adiabatic regime. 

The phonon PDF can be easily derived within LF ap- 
proach. Being the local electron densities parametric 
variables in the anti-adiabatic limit the phonon PDF can 
be written as 



P(X) =^ W/ P (*- X l 



(38) 



where wi is the probability of having an occupancy m, 
Xi the relative displacements and Po(X) the ground state 
PDF of an harmonic oscillator. Pq(X — Xi) is then the 
conditional probability of having a displacement X given 
a definite occupation m. In the spinless case n; = 0, 1 
with equal probability giving 

P{X) = I- (P (X - X ) + P (X + X )) (39) 



where Xo = ia/y/2, A definite polarization can be asso- 
ciated to the ground state if the PDF becomes bimodal. 
By requiring dP{X)/dX\x=o > 0, which simply means 
that the X = turns into a local minimum from the 
maximum it is in weak coupling, we get the usual anti- 
adiabatic condition for the existence of a polaronic state, 
i.e., a 2 > 1 (see fig. EJ. 

In the spinful case ni = 0, 1, 2 

P{X) = n d (P Q (X-2X )+P Q (X + 2X Q )) + {l-2n d )P (X) 

(40) 

where n d — (n^n^) is the site double occupancy. It is 
worth noting that in the insulating state n d — 1/2, and 
the zero-displacement PDF associated to singly occupied 
sites is depleted. This is an example of the strong cou- 
pling dependence of the phonon properties on electronic 
state. In contrast with the spinless case the existence of 
a definite polarization is now associated to a bipolaronic 
state. The condition under which (|40J) becomes bimodal, 



exp(-2a 2 )(4a 2 - 1) > 



l-2n d 
2n d 



(41) 



An estimate for the bipolaronic transition can be ob- 
tained by taking n d — 1/2, which gives a 2 > 1/4. The 
presence of a fraction of singly occupied states increase 
the critical value of a needed for a bipolaronic states to 
be formed (see Fig. We notice that the spinful PDF 
for rid = 1/2 maps onto the spinless one after the usual 
rescaling (A/2 — > A and 27 — > 7). 

Finally let us compute the phonon propagator in the 
anti-adiabatic regime. Using H30|l we obtain 



D{t)=D (t)+4a 2 sE d {t) 



(42) 



where Dq (t) is the correlation function of the harmonic 
oscillator of frequency loq and 



H d (t) = -i/T^K(t)-i)(n 



(0)4) 



(43) 



is the density-fluctuation correlation function^ 3 .. Let us 
first examine a limit case in which both the spinful and 
spinless fermions are insulating, i.e., the atomic (zero hy- 
bridization) case D = 0. In this limit the density is a 
constant of motion and the Fourier transform of 14211 is 
we obtain in the frequency domain: 



D{u) = D q (uj) - i^a 2 - A 
D(lo) = D q (uj) - iAa 2 (n^ni) 



1 ( 1 



1 



4 \ ui + iS to — iS 
1 1 



u> + iS ui — iS 



(44) 
(45) 



respectively in the spinless and in the spinful case. From 
Eqs. and (|43|) we see that the pole associated to the 
density fluctuation occurs at zero frequency. In other 
words the freezing of the density fluctuations which oc- 
curs at strong coupling imply a vanishing of a low fre- 
quency contribution in the phonon DOS. Of course all the 



13 



phonon spectral weight remains at frequency ±wo- Upon 
including the hopping term, the situation changes dras- 
tically depending on the occurrence of a MIT. If a MIT 
takes place, the vanishing of the zero-frequency spectral 
weight is attained at finite coupling (Xmit) on the other 
hand in the spinless case this occurs only asymptotically. 
To corroborate such a prediction which explains the be- 
havior shown in Fig. [21 we perform an approximate cal- 
culation of the correlation function l|42|) . which gives for 
the spinless case and in Matsubara frequencies, 2^ reads 



S d (iw n ) = - G p (iu>„ 



iLu n ).G p {iuj m ) (46) 



If we consider the motion of the spinful electron within 
the CPA approximation as a motion in a random bimodal 
potential^, then the same expression holds also in the 
spinful case, and the differences are only provided by the 
form of G p (iuj m ) ■ 

Introducing the spectral representation of G pi perform- 
ing the sums and the analytical continuation to real fre- 
quencies we obtain 
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FIG. 14: (color online) Phonon DOS in the spinless (upper 
panel) and spinful (lower panel) cases in the anti-adiabatic 
regime (7 = 2.0 spinless, 7 = 4.0 spinful). Solid curves are 
DMFT results, dashed lines the LF for the approximation 
for the charge fluctuation contribution in Eq. (11 12 i . The 
values of A are in the spinless (spinful case) from left to right: 
A = 0.4 (0.2), A = 1.2 (0.6),A = 2.0 (1.0),A = 5.2 (2.6). 



dvp p (v)p p {v - lu) 



(47) 



where p p (w) = —ImG p (uj + iS)/n. Notice that if we 
consider the atomic limit t = we recover l|44H . but not 
1)450. This is a drawback of the CPA approximation for 
the density correlations. We adopt however Eq. H47|) 
as a CPA prescription which qualitatively predicts the 
disappearance of the phonon spectral weight at the MIT 
in the spinful case as we shall see explicitly below. 



B. Discussion of the results in the anti-adiabatic 
regime 

In Fig. ^|we compare the phonon DOS from DMFT- 
ED with the LFCPA described above. A strong reso- 
nance of phononic origin is always present around too, 
together with a low-energy broad structure due to elec- 
tron density fluctuations predicted by HLFA ||32J. 

As in the adiabatic case, but here even more clearly, 
we see that the MIT implies the disappearance of low- 
frequency weight in the phonon DOS. The disappearance 
takes place exactly at the MIT, i.e., for a much smaller 
coupling than bipolaron formation. In a sense, the effect 
of the electrons on the low-energy part of the phonon 
spectrum is always similar to the antiadiabatic regime. 
A noticeable difference between the present antiadiabatic 
case and the adiabatic regime is instead present in the 
high energy part of the phonon spectrum, which in this 
case is basically unperturbed in the present antiadiabatic 
regime. We notice incidentally that, unlike the present 
case, phonon softening has actually been found in spin- 
less Holstein model in one dimensions'*'. There how- 
ever a MIT is induced by Peierls dimerization and the 



phonon softening is a precursor of the MIT just like in 
our adiabatic case where A s < Xmit- 

In our previous publication*^ we have found a soften- 
ing which is not complete even in the spinful by esti- 
mating the renormalized phonon frequency from the 
uj n = limit of the phonon self-energy in the imaginary 
frequency (f2 2 = ujq + 2u} H(iuj n — 0))). We observe here 
that this quantity is a good estimate for the phonon fre- 
quency only when the majority of the phonon spectral 
weight is concentrated around some frequency, and the 
spectrum is not divided in several features. A straight- 
forward calculation gives 



1 [°° 1 



(48) 



which implies that all the different features of the phonon 
propagator contribute to f2, therefore leading to an in- 
complete softening also in the case in which the lowest 
phonon pole actually completely softens. 

In Fig. E| the electron DOS is compared with the 
results of HLFA. The different behavior of spinless vs 
spinful cases is marked here clearly by the presence of 
the MIT in the former at a value of the coupling which is 
much less that that of the (bi)polaron crossover. HLFA 
correctly catches the gross behavior of the DOS. Notice 
that at the strong coupling the CPA employed to obtain 
the lower diagrams accurately reproduces the position of 
the side bands. 

The different behavior of spinless vs spinful system 
can be easily understood in terms of strong coupling 
anti-adiabatic perturbation theory for the original lattice 
problem 37 -. Here we sketch an alternative approach based 
on the equivalent HIM J2J) after the LF transformation 
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FIG. 15: (color online) DOS in the spinless (upper panel) 
and spinful (lower panel) cases in the anti-adiabatic regime 
(7 = 2.0 spinless, 7 = 4.0 spinful). Solid curves are DMFT 
results, dashed lines the LF approximation. The values of A 
are the same of Fig. 1141 



(|32[l In the spinful case at second order in hybridization 
V k the following Kondo hamiltonianS can be obtained^: 



H = H' 

+ 4 



II 



VkPfPc(k) 



J2 v MPc(k)+Pjpt(k)) (49) 



where H 1 contains all the terms of 132(1 which are not 
proportional to V k , V 2 = £ fc \V k \ 2 and \v k \ 2 = \V k \ 2 /V 2 . 
p are pseudo-spin operators corresponding to bath (p c ) 
or to impurity (/?/) states, and J's are anisotropic Kondo 
couplings (see also Eq. (8) of RefHoT) 



J11 



8V 2 
D 



E(±) 



Q 



2m 



m!(m7 + A/2) 



(50) 



where the + (— ) sign is taken for the || (_L) coupling. 
In the spinless case the processes leading to J± do not 
exist while the remaining Ju is solely associated to charge 
fluctuations i.e. it will correspond to a /-charge c-charge 
interaction in the effective Hamiltonian 



H = H' + J l] J2vkP Z fP z c (k) 
fc 



(51) 



where now P% = y f — 1/2. No Kondo effect is expected 
for this Hamiltonian, while in the spinful case the model 
will display a Kondo effect associated to the pseudo-spin 
fluctuations. A strong coupling estimates of Ju gives 



J 



8V 2 
DX 



(52) 



and J±/J\\ oc cxp(— 2a 2 ) which means an exponential 
suppression of the superconductivity versus charge cor- 
relation at strong coupling due to retardation effects^. 



In the anti-adiabatic limit 7 — > 00 we recover an isotropic 
Kondo Hamiltonian J± = Ju . 

It is worth to note that the range of validity of these 
expansions is somewhat larger than the anti-adiabatic 
regime. Indeed the minimum energy scale of Hamilto- 
nian H' should be larger that the largest Kondo coupling 
just evaluated i.e. Ju at strong coupling. Therefore sup- 
posing ojq to be smaller that the bandwidth of the "bath" 
electrons and of polaronic energy we get 7 > 4V 2 /D 2 X. 



VI. CONCLUSIONS 

We have presented a thorough study of the Holstein 
model comparing the spinful and the spinless Holstein 
model at half-filling with the main purpose of identify- 
ing in the most unambiguous way the polaron crossover, 
and disentangling the tendency to form polarons from 
the metal-insulator transition associated to bipolaronic 
binding. The role of the phonon frequency, determin- 
ing whether the system is adiabatic or antiadiabatic 
(slow/fast phonon dynamics) has also been addressed 
with care. We computed various observables within Dy- 
namical Mean-Field Theory, solving the associated im- 
purity model by means of Exact Diagonalization. By 
comparing these numerically exact results with approxi- 
mate analytical schemes suitable for the different limiting 
regimes we have identified different mechanisms for the 
appearance of polaronic features in the different regimes. 

In the adiabatic regime a Born-Oppcnhcimcr scheme 
has been devised to deal with both small and strong cou- 
pling regime. In the anti-adiabatic case a more standard 
Lang-Firsov approach for the impurity model has been 
applied. 

While in the adiabatic case the polaron crossover and 
the bipolaronic MIT occur for similar values of the cou- 
pling, in the anti-adiabatic case pairing occurs for smaller 
coupling than polarization. Therefore there is a large 
region of the phase diagram (whose size increases with 
the phonon frequency) in which pairs are formed, but 
no lattice polarization is associated. At the bipolaronic 
MIT the vanishing of quasi particle spectral weight at the 
Fermi levels is accompanied by a softening of the phonon 
mode associated to low-energy charge fluctuations. The 
absence of such a complete softening in the spinless case 
is a further signature of the absence of insulating behav- 
ior. 

Within DMFT the Mott-Hubbard transition is associ- 
ated to the Kondo effect of the Anderson impurity model. 
Analogously the pairing transition of attractive models 
is related to the Kondo effect once a pseudospin whose 
components are the s-wave superconductivity and charge 
ordering are introduced. In the Holstein model we have 
a mapping onto an anisotropic pseudospin Kondo model, 
in which pairing correlations are greatly reduced by re- 
tardation effects but nonetheless survives due to phonon 
quantum fluctuations. At the MIT the Kondo peak as- 
sociated to quasiparticle properties shrinks to zero and 
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phononic low energy features disappears. 

An interesting feature of the polaron crossover in both 
spinful and spinless case is the failure of the BO ap- 
proximation in the crossover region. At the crossover 
phonon and electron state becomes entangled leading to 
a breaking of the BO approximation even in the adia- 
batic regime. We plan to extend the Born-Oppenheimer 
scheme here introduced to the more involved situation in 
which the electrons also feel a Hubbard-like repulsion, a 
subject which has been the subject of different numerical 
analyses, using ED^i^Sil^, and NRG techniques^. 

The authors acknowledge useful discussions with E. 
Cappelluti and C. Castellani. This work was supported 
by INFM PRA-Umbra and MIUR-Cofin 2003 matching 
funds programs. 



APPENDIX A: PERTURBATION THEORY OF 
THE TSPM 

1. Spinless Case 

The spinless case is obtained by taking s = 1 in 
(|14fl . We will treat the hybridization term of (|14|) as 
a perturbation of the free Hamiltonian defined by Hq — 
-2e (/+/ - ~) a z - Aa x + J2k e kc\c k . We take the la- 
bel k to assume symmetric values with respect to zero. 
We impose the particle-hole symmetry at half-filling by 
means of the following conditions: 

V k = -V. k 

£fe = -e-fc. 

If conditions ilAlfl are obeyed, Eq. (|14f) is invariant for 
the following transformations: 

^± 



(Al) 



--k 



where a = ±1. The free propagator of the impurity is 
defined as: 

Go(t) _ _ Me--W W ) , m 



Using (|A3|) we get 



1 



, cosh(2rZ - pi) e 2 A 2 



1 



-r) 



cosh(2rZ + /3l) 



A 2 



2" v ' ' \ cosh/3/ I 2 I 2 
which, after a Fourier transform, becomes 



IT > (A5) 



G o (i0 n ) = i^G OQ (i0 s ), 



where we have introduced Go a given by 
„ ,. s 1 A 2 1 A 2 1 



2 \l 2 i0 n + 2la I 2 i0 n 



(A6) 



(A7) 



and 0„ = (2n + l)n/p. 

The hamiltonian (|14H is quadratic in the fermions but 
the spin operators generate an effective interaction. For 
this reason we can not expect a finite set of equations of 
motion if the band contains "many" electrons. 

We introduce the quantity : 

G ff+ (r) = -|tr(e-^T r /(r)/+(0)) , (A8) 
we find: 

G ff+ (r) = -5(r) + 2eG^ ff+ {r) - $> Cfe/+ (r), (A9) 

fc 

where adopting the same notations as for Eq. l|A8jl . 

G<j z ff+{ T ) is defined as 



where (/~ = /, cr = c). It is easily verified that the 
previous transformations are generated by the unitary 
operator U = exp[(6/ - f + b+ + J2 k ( d -kCk - d± k ) + 
i&x)'K/2] (b and d k are the new particles). Of course the 
total particle number N = /+/ + J2k c t°k 1S conserved. 

Hq is easily studied. This Hamiltonian is diagonal in 
the fermions. It is immediate to treat the band electrons 
term. The spectrum and the eigenvectors of Hq are given 
in Eqs (|17fl and (|18|) . Now we introduce 



1 



A 

e — la 



y/2l(l - ae 

1 fe + la 

A 



y/2l(l + ae) 
The following properties hold: 



(A2) 



A 



(A3) 



G. zfS+ {r) = --tr ( e -^T T ^(r)/(r)/+(0)) . 

G Ck f+ is defined similarly. The equations of motion for 
the r.h.s. terms of Eq. (|A9|) are 



G Ckf+ (r) = -e k G Chf+ (r) -V k G ff+ (r); 
G azff+ (t) - iAG auff+ (r) + 2eG ff+ (r)- 
-Y^V k G^ Ckf+ (r). 



The first equation does not give rise to any new quantity, 
but the second involves G rTy ff + and G azCk f+. Therefore, 
we should derive the equations of motion for these terms, 
and proceed until the iterations do not generate any new 
quantity. Unfortunately this procedure does not seem 
to be convergent. Terms like G ayCk f+, G ITx ff+, etc. are 
generated and the equation for G ayCk f+ introduces the 
correlator 



a=±l 
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G 



z^e fc /+(T) = --He-^T T 



-2e[r(r)f(T)--)a x (r)c k (r)r(0) 



The equation of motion for the latter involves, among to resort to perturbation theory. 



others, 



-tr { e-' 3H T T 2e ]T V p c+{r)f{r)a x (r) Cfe (r)/+(0) 



We aim to compute the interacting propagator: 
tr(e-^T r /(r)/+(0)) 



G(t) = - 



tr (e-P H ) 



It is evident how a hierarchy of correlators with increas- We resort to second order perturbation theory in the Vk 

ing complexity is thus generated by this process. The coefficients. If Zq = tr(e~P H °), V is the hybridization 

impossibility of further pursuing this approach forces us and, as usual, V(t) = e TH °Ve~ TH ° , we have: 

I 



1 







>z -tr {e^ H «T T f(T) J dndriVinW^f+iO) 



/ connected 

Taking the Fourier transform and using again Eq. I|A3(1 . we obtain the second-order perturbative result 



1 



a 

r \ 



v 2 

E k 



I0n - £fc 



1 



eA 

I 2 \i0n + ^la 10 



-y 



atanh(/3Z) 



, iPn - Cfe i0n - iPn + 21 a ' 



r 



(A10) 



where 0„ and p n are fermionic momenta. The interpre- 
tation of the first term of Eq. (|A10(1 is straightforward 
since in can be written as 



G 



A^i^n) = ^yG Oa (i0 n )T, ihl (i0 n )G O a(i0n), (All) 



Eq. (|A6|) , can be interpreted as the free propagator of the 
impurity for a given spin a. As a first approximation for 
the interacting propagator we can sum the Dyson series 
for the j4-terms, obtaining the BOCPA approximation 
Plfl discussed in Sec. |IV| 



where we have introduced the hybridization self energy 



Sibr(«0n) = J2 



k i0 n 



Goct, defined in the last line of 



The second term (henceforth called £?-term) of Eq. 

Engl, 



Gg } ( ? 0„) = ^ 



eA 



1 



i0 n + 2k 



1 

i0n 



E 



v k 2 



atanh(/3Z) 



£fc i0 n - ip n + 2la ' 



(A12) 



is surprising for the following reasons: 

i) it does not contain the free propagator; 

ii) the two external legs, namely the two factors 
carry a fermionic frequency i0 n but 



L.) 



the sum of the residues of the poles is zero so this object 



qtanh(ffl) 



carries a bosonic fre- 



does not carry charge; 
Hi) the internal line . 

7 i0 n — ip n +2la 

quency and, when summed over a, again has residues 
adding up to zero. 

The last two properties suggest us to interpret the 
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propagators of B-type Gjp as spin correlators. 

To put this interpretation on more solid ground, we 
consider the spin correlator defined as: 



S'*(r) 



tr„ (e-^°T T a z (rK(0)) 



tr Q (e~f 3H °) 
which is explicitely 

E a K\e- fm °T T a z (r)a x (0)K) 



S z a x {r) 



(A13) 



(A14) 



The imaginary time ordering T T is defined taking 
into account the canonical anticommutation relation: 
Wz,(?x} = 0: 

T t o z (t)o x (t') = 6{T-T r )a z {T)a x {t')- 
- 9(t' -T)a x {t')a z {T). 

The r — > t' ordering limit can be defined arbitrarily ac- 
cording to convenience. It is easy to find that S^ x (t) is 
antiperiodic in r: if r < but r + [3 > then 

Of course only odd frequencies (0 n = |j(2n + 1)) con- 
tribute to the Fourier series of S^ x (t): 



S Z a X {i0n) 



eA / 1 

I 2 \i0 n + 2la 



1 

i0r, 



(A15) 



Therefore, the spin-correlator reproduces the first factor 
in Eq. I|A12|) . It is worth noticing that if we consider the 
correlator: 



S x a z (r) 



tr a (e-f* H °T T (j x (T)(T z (0)) 
tr Q (e-^o) 



we obtain the same result (|A15|) (the equal-time ordering 
of a x and a z is irrelevant). 

We still have to explain the last factor of Eq. I|A12(I . 
Now the product a z <j x is proportional to a y , so it is nat- 
ural to consider the correlator: 



S*»(r) 



tr Q (e-P H <>T T a y (T)<T y (0)) 



tiv, 



-PH \ 



(A16) 



It is not obvious how to define a time-ordering operator 
in this case. Since = 1, a fermionic time ordering does 
not seem appropriate (/ 2 = for a fermion). For this 
reason we define a time ordering taking into account the 



canonical commutation relations: [cr y ,a y 



0: 



T T a y (r)a y (r') = e(T~T')a y (r)<j y (t')+0(r' -r)a y (t')a y (r). 

It can be observed that 5™ is actually independent on 
a, so we will omit this index. For the Fourier transform 
we find: 



S VV (i0 n ) 



Ea tanh(/3Z) 
i0 n + 2la ' 



(A17) 



with bosonic frequencies (0„ = ^S-n). This is exactly the 
missing piece of information to complete the interpreta- 
tion of Eq. ijSnjl. 

Two useful results are 

Qz*U* ^ A 2 \ - atanh(ffl) 
S ^)=—pl^ l0n + 2 la 



S XX (i0 n ) 



e 2 ^ a tanh(/3Z) 
I 2 ^ i0 n + 2la ' 



where S xx and S zz are defined analogously to S vv and 
the frequency n is again bosonic. 

In order to define the diagrammatic expansion used in 
Fig. IS] we need four types of lines (propagators) and two 
vertices. The lines are: 

e 2 1 A 2 1 

impurity: (/+/)„= w - + 



I 2 i0 n + 2la I 2 i0„ 



bath: 



1 



(Cfc c k ) = 



I0n - £fc 



/-spmon: {a z a x ) = — 



I 2 V i0n + 21<X I0 r . 



b— spinon: (a y a y ) 
The vertices are: 



atanh(/3Z) 
i0 n + 2la 



-H 



where each bold triangle denotes the hybridization con- 
stant Vfc. 

Let us detail the Feynman rules (derived from the 2nd 
order result). We start from the conservation laws. Each 
line carries a spin and a fermion number. The fermion 
number nj is defined to be 1 for Ck , / and f-spinon lines 
and for the b-spinon line. The spin number n a is 1 for 
the spinon lines and equals for the / and Ck lines. 

(a) Fermion and spin numbers are separately conserved 
in every of Eqs. I|44I45|) . allowed process. For the pur- 
pose of computing G(i0 n ), the band electrons and the 
b-spinons act like virtual particles. While this observa- 
tion was relatively expected, a much more striking result 
of Eq. i|A10|) is that the f-spinon lines can also be real: 

(b) The diagrams contributing to G(i0 n ) are obtained 
drawing self energy diagrams with external lines 



According to this rule the f-spinon line should con- 
tribute to zero order, whereas it does not. The reason 
is that a sum over the a index is actually implicit and, 
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according to Eq. I|A15(1 . this sum vanishes. As a conse- 
quence of (a), self-energy diagrams with external legs of 
different types are forbidden. 

(c) All the lines contributing to a given diagram have 
the same a label (the a label should not be confused 
with the n a number). We have to sum over a, with a 
coefficient 1/2 if the two external legs are (free) impurity 
propagators, and 1 if they are f-spinons. (One may try to 
justify this rule taking into account the sum over the a 
index: the impurity selects a value for this index, depend- 
ing on the time ordering, while this is not the case for 
the spin operators, according to Eqs. (|A13I) and IjAlSf) ). 
These rules lead to diagrammatic representation of the 
perturbation expansion shown in Fig. |H1 



2. Spinful case 

The eigenstates and eigenvalues for Hq = 
~ € Y, a =± (Ufa ~ \) a* ~ Acr^ are easily found. 
Denoting with the singly occupied states and with 
\u) the doubly occupied or empty states we get eigen- 
vectors and eigenvalues reported in Eqs.. (|20|l . (|21|l and 



Using the same notation of Eq. (|A2|) we get 

A + aa(l — ae) 



2y/l{l - ae) 
A — aa(l + ae) 



-a- 



where a, a = ±1. The scalar products (|A18(1 do not 
depend on a' . The completeness relations are: 



l"S)(«Sl = i; 

a=±l 

]T \<)(<\ = i. 

a=±l 



Let us define the following propagators: 



(A19) 



Gla 



I- A 



l + A 



1 

21 \i0 n + a(l + A) ' i0 n + a(l-A) 
1 / l-A l + A 



21 \i0 n + a(l + A) i0 n + a(-l + A) 



(A20) 



The free propagator is given by: 



G o (i0n) = V (cosh(/?Z)Gi Q + cosh(/3A)G 2a ) . 
Zi + — ' 

(A21) 

where Z + = 2cosh(/3Z) + 2cosh(/3A). The second-order 
perturbation theory in Vk displays a structure similar to 
the spinless case. After a straightforward calculation we 
find, including zero order term: 



G [2) {i0 n ) = — ^] (cosh(pi)[G la (i0 n ) + G la (i0 n )Y; ihl (i0 n )G la (i0 n )] + 

+ COsh(/?A)[G 2a (i0„) + G 2 a(i0™)Sibr(i0n)G2 Q (i0„)] 



1 1 ^fe e / a silL 

Z+ a P iPn _ £fe 12 ~ P 



asmh(f]l) 



"(- 



-p n ) + 2la y i0 n + a(l + A) i0„ + a (I - A)' 
asinh(/3A) , 1 1 



-p n ) + 2Aa y i0 n + a(l + A) i0 n + a(—l + A) 
I 



n 



(A22) 



As for Eq. (|A10(I . the two terms are of diffent type. 
In particular the free propagator does not appear in the 
second term, the external lines do not carry charge, and 



inside the loop there is a chargeless bosonic line. As for 
the spinless case, let us consider the Dyson series for the 
first type terms (coherent potential approximation): 



G c {i0 n ) 



1 ^ / cosh(/3Z) cosh(/3A) 

Z+ \Giaii0n)- 1 ~ S ib r(i0n) G 2Q (i0, i ) _1 - Sibr(«0n) 

I 



(A23) 



Since I > A, in the zero temperature limit (3 — > oo, we get for this approximation of the propagator: 

1 



Gc{i0n) = \ Y1 



2 ^ Gi a (i0 n ) 1 - Eibr(«0 n ) 



(A24) 
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which is the BOCPA in the spinful case. 
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